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Primer vector theory u used in analyzing a set of linear, relativ e n o tion equa- 
tions - the Clobessy-ULltshire (C/M) equations - to determine the criteria aid 
necessary conditions for an optimal, l-lnpulse trajectory. Since the state 
vector far these equations is defined in terns of a linear system of ordinary 
differential equations, all fundamental relations defining the solution of the 
state and oostate equations, and the necessary conditions Tor optimality, can be 
expressed in terns of elementary functions. The analysis develops the analytical 
criteria for improving a solution fay (1) moving any dependent or independent 
variable in the initial and/or final orbit, and (2) adding intermediate impulses. 
If these criteria are violated, the theory establishes a sufficient number of 
analytical equations. The subsequent satisfaction of these equations will 
result in the optics! position vectors and tines of an ■ -Impulse trajectory. 

The solution is examined for the specific boundary conditions of (1) fixed-end 
conditions, two-impulse, and time-open transfer; (2) an orbit- to-orbit transfer; 
and (3) a generalized rendezvous problem. A sequence of rendezvous problems is 
solved to illustrate the analysis and the computational procedure. 


2.0 DtTBODOCTIOW 

In the early 1960's, the solution of an optimal trajectory, based on the assump- 
tion that the thrusting acceleration is replaced by an impulse, received consid- 
erable attention as a method of performing mission analysis studies. Lawden 
(ref. 1) developed the necessary conditions for an optimal impulsive trajectory. 

He examined the limiting conditions on an optimal finite thrust solution where- 
in the thrust magnitude was unconstrained but bounded between a maximum and mini- 
mum value. These results are known as Lawden' s necessary conditions for an opti- 
mal impulsive trajectory. They specify the conditions that must be satisfied by 
the primer vector and its derivative on a candidate impulsive trajectory. 

Lion and Handelsman (ref. 2) later established the criteria and necessary opti- 
mal conditions for a fixed-time, impulsive trajectory whereby a reference 
trajectory oould be improved (i.e., the cost function or the sum of the magnitudes 
of the applied impulses could be decreased). 

Jezewski and Rozendaal (ref. 3) combined the Lion and Handelsman (ref. 2) re- 
sults with a oonjugate gradient iterator (ref. M) to produce an algorithm that 
generates optimal impulsive trajectories in an efficient and rapid manner. 

In the middle 1970's, Jezewski (ref. 5) extended this analysis to generate 
a general differential cost function for the two-body problem. This function 
defines the gradient structure and cost function for (1) any set of boundary 
conditions when the applicable constraints are specified, and (2) equality 
and inequality constraints on both the state and control variables. By this 
means, completely general, two-body, N-impulse, optimal trajectories can 
be generated for any set of constraints that can be expressed mathematically. 

This study applies primer vector theory to a set of relative motion equations 
- the C/W equations - to determine the necessary conditions for an optimal, 
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■-impulse trajectory. In previous studies, Prussing (ref. 6) applied primer 
vector theory to the problem of obtaining Multiple impulse, fixed- time rendez- 
vous solutions between coplanar circular orbits. Jones (ref. 7) used an elanent 
formulation to determine a general optimal set of lnpulses far every rendezvous, 
as loc^s as the linear assumptions were not violated. These and other similar 
studies construct priner vector rendezvous solutions that satisfy the necessary 
conditions by choosing and iterating directly on the constants of the solution. 

In this study, completely general, optimal, M- Impulse transfers between given 
boundary conditions are obtained for the linear system of equations. The gen- 
eral rendezvous problem is just one subset of these solutions. The technique 
for achieving a solution differs from those previously mentioned and is similar 
to the method outlined in reference 5. A general differential cost function is 
developed in terms of a set of all possible independent parameters. The spe- 
cific problem to be solved is selected from this set; the coefficients of the 
specific set of parameters are required to be zero to achieve an optimal solu- 
tion. This general method has proved to be highly successful in the nonlinear 
problem (ref. 5). ?' 

In the primer vector theory, the solrMon of the state and costate equations is 
used along with a scalar function known as the adjoint equation to develop the 
criteria and necessary conditions for a candidate trajectory to be optimal. 

Since the C/M differential equations are a linear system of ordinary differen- 
tial equations, they can be analytically integrated in terms of elementary 
functions. Also, the costate differential equations (the costate vector is 
canonically conjugate to the state vector) are also linear and can be integrated 
in terms of elementary functions. Hence , all fundamental relations defining the 
state, costate, and necessary conditions for optimality are known analytically 
in terms of elementary functions . 

Section 3.0 will proceed to define the differential equations for the state and 
costate vectors in a convenient form and to integrate these equations in terms 
of elementary functions. A scalar function known as the adjoint equation will 
be developed from these differential equations, which will be used in the subse- 
quent development of the necessary conditions for optimality. 

In section 4.0, a cost function is defined and, after taking its variation and 
using the adjoint equation, a general differential cost function is developed 
for the C/W equations. This function defines the necessary conditions to be 
satisfied by the C/W equations for an N- impulse trajectory for any changes in 
(1) the initial »nd final orbits, (2) the initial and final times, and (3) the 
times and position vectors of any intermediate impulses. The criteria for an in- 
termediate impulse is developed in section 4.1, and subsequent analysis is 
performed to determine the position vector for this impulse. 

Finally, the general differential cost function is examined for the specific 
boundary conditions of (1) fixed-end conditions, two- impulse, time-open trans- 
fer; (2) orbit-to-orbit transfer; and (3) a generalized rendezvous problem. 

A number of example rendezvous problems are solved to illustrate the analysis. 
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3.0 APPLICABLE EQOATIQRS 


3.1 STATE VECTOR EQUATIONS 

The development of the linear, relative-notion differential equations can be 
obtained fron reference 8. A particular fora of these equations, know as 
Hill's equations (ref. 9) or the Clohessy-Viltshire (C/V) equations, are: 


x - 2wfr = 0 

y + 2wi - 3 u> 2 y = 0 

z ♦ s 0 (1) 


where w is the constant angular velocity of a coordinate system, ass used to be 
located in a circular orbit of a given radius magnitude. These equations, which 
describe the notion of a particle with respect to this coordinate systen, can be 
analytically integrated such that the solution is described in terns of the 
boundary conditions on the trajectory and the solution tine. 

Defining a state vector as 


S T = (R t ,V t ) 


( 2 ) 


where 


R t = (x,y,z) , V = R 


and the superscript T refers to the transpose. The solution to equation (1) 
can be expressed as 


S(T) = A(t,0) S(0) 


(3) 


where the relative time t is the difference in the absolute times t (for ex- 
ample, T i * t 2 - t-| ) and the matrix A (a function only of w and t) maps the 
state at one time into another time. If ayj(i,j = 1,2,..., 6) are the elements 
of this matrix, then the only nonzero terms are 


80FH5 


*11 * 1 

*22 

*12 = 6(urr - s) 

*24 

a^jj = (4s - 3ut)/<i> 

*25 

a^ s 2(1 - c)/(ii 



*52 

ai |2 - 6q)( 1 - c) 

*54 

aj|i) = 1,r - 3 

*55 

ajjg = 2s 



4 - 3c 

*33 

s C 

2(c - 1)/« 

*36 

= S/(0 

3/(0 




*63 

s -d)S 

3«s 

*66 

S C 

-2s 




s C 


where 3 and c are the sin(arr) and cos(urr), respectively. If we write 
equation (1) in first-order form using equation (2), the C/W differential 
equations can be expressed as 


S = F(R,V) 


(4) 


where the vector F is defined a3t 


F^ = (v x ,Vy,v z ,2 a) Vy,3a} 2 y — 2( 1 )V x ,-q)2z) 


(5) 


Equation (4) will be used in sections 3.2 and 3-3 to develop the costate 
equations and the adjoint equation, both of which will be used to subsequently 
develop the necessary conditions for an optimal N-impulse C/W trajectory. 


3.2 COSTATE VECTOR EQUATIONS 

The costate equations, or the Lagrange multiplier equations, are a system 
of equations adjoint to the state equations. In differential form, they 
are defined as 


X 



( 6 ) 


writo.-c *-he vector F is given in equation (5). 
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Developing these differential equations, we obtain 


^1=0 
^2 s - 3^2X5 

*3 = ^Xg 


^l( = ♦ 2^5 

^5 s -^2 - 2*^4 

*6 = -X 3 


Since these equations are linear, they can be immediately integrated to obtain 


Q = 


P = 


-<* 1 / 3 U 

o^c + O3S + 2(014 O^t) 
10(053 - Ogc) 


-2 

— (o^c + 03s) - (04 + ajt)/u 

(o^s - a 3 c - 20^0)) /310 
a^c + 05s 


( 7 a) 


( 7 b) 


where the vectors Q and P (primer vector) are defined by the relation 


X T = (Q T ,PT) 


( 8 ) 


and the constant coefficients c^(i = 1,2,..., 6) are boundary conditions 
computed as follows. From reference ( 5 )» the primer vector at the time of 
an impulse is defined as a unit vector in the direction of the impulse, or 



( 9 ) 


If we designate Pi T = P T (ti) = (11,12,13) and similarly, 

P2 T s P T (t2) = (mi, 1112,1113) , then the in-plane elements of the coefficient 
vector o in equation (7) are determined as 
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’ tt 1 


» « 
ll 

a 2 

= B" 1 

l 2 

a 3 


“1 



m 2 

■ m 


where the only nonzero elements bjj(i,j = 1,2,. ..,4) of the matrix B 
(determined from eq. (7b)) are 


b 12 

= -2/3w 

b 21 

= bi 2 /a) 

bm 

= — 1/U) 

b23 

= -1/3a) 

b 3 i 

= -T/u) 

b 4l 

= b 21 

b 32 

= -2C/3 a> 

b 42 

= -b 23 s 

b 33 

= -2s/3w 

b4 3 

= b 23 c 


b 3 n = bin 

Since the differential equations for the state vector are uncoupled in the 
z-direction, the out-of-plane elements of the coefficient vector a are 
determined as 


a 5 s x 3 

05 = (m3 - l 3 c)/s 


( 10 b) 


Note that the out-of-plane solution is not valid for u)T r nir, n s 0,1,... 

Knowing the value of the vector a , the primer vector P and its companion 
vector Q can be determined analytically from equation (7) for any time t. 

The costate differential equations could also have been obtained from a 
Hamiltonian approaoh. If we define a Hamiltonian H as 


H S P T V + qTr 


( 11 ) 


then the state and costate C/W differential equations can be obtained from 
the canonical form 


6 
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RT 


3H 
* 3Q 


VT s 


3H 

3P 


and 


QT s 


3H 

3« 





In later sections, we will have occasions to use not only the concept of a 
Hamiltonian, but also its mathematical definition (eq. (11)). 


3.3 ADJOINT EQUATION 


% 

We need one further mathematical relationship before we can prooeed with the de- 
velopment of the criteria for improving a C/W trajectory and the neoessary condi- 
tions for an optimal, N-impulse solution. A functional relationship is required 
between the oostate vector and the perturbations, in the state. 


From equation (l|), the differential equation for the variation in the state 
vector can be expressed as 


N 


6S * (jlJ* S (12) 


where £ represents the contemporaneous variational operator. Note that 
because the vector F is linear in the state veotor, equation (12) represents 
the exaot variational differential equation and not a truncated expansion. 
Premultiplying equation (12) by \ T and equation (6) by £S T and adding, 
we obtain 


X t 6S + fiS T X = X T 6S - £S T (— ) X 
\3S/ \3 S/ 


The right hand side of this equation oan be easily verified as zero. The 
left-hand side of this equation oan be expressed as the exact differential, 


d 

dt 


(X t 6S) * 0 
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which implies that on any C/W arc, the scalar product of the coatate vector and 
the variational state vector is a constant, or 


X t 6S = constant 


(13) 


Equation (13) will be referred 


to as the adjoint equation for the C/W 



tions . 


11.0 COST FUNCTION AND THE NECESSARY CONDITIONS FOR OPTIMALITY 

The cost function, J, is defined as the sum of the magnitudes of the applied 
impulse vectors. Let us assume that the solution consists of three impulses 
with the following notation: at the time t-|, the impulse vector is AVj; at 
the time t 2 , the impulse vector is AV 2 ; and at an intermediate time tgj, .0 
(t^ < t m < t 2 ), the impulse vector is AV m . A three- impulse solution has been 
assumed for the analysis because the results can be (1) easily reduced to a £ 
two- impulse solution by the removal of the intermediate impulse, and (2) readiljr 
extended to more than three impulses by additional intermediate impulses. This 
three- impulse cost funotion is 


J = |A»l| ♦ |AV m | ♦ |AV 2 | 


(14) 


where AV = V + - V" and the superscribt minus and plus refer to the evaluation 
immediately before and after an event, respectively. To determine the necessary 
conditions for optimality, we shall examine a particular perturbation - a con- 
temporaneous variation - of the cost function. Taking the variation of equation 
(14; t we obtain 


6J 


AVi t 


6 (AVt) 


av b t 

lAVmi 


5 (AV m ) 



(av 2 ) 


Using equation (9), evaluated at the times t-j, t^, and t 2 
and taking the variation of the individual impulses, we have 


in this equation, 


6J * P^fiV^ - 6 V,-) «■ P m T (6V n + - 6V) ♦ P 2 T (6V 2 + - «V 2 ") (15) 


The terms P^Vi* and P 2 ^6V 2 " can be eliminated from this equation by 
evaluating the adjoint equation (eq. (13)) on the two arcs separating the im- 
pulse : 


8 
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(Pi, «v^) * «^-, «v> ♦ <p*“, «V> - <Qi. «*i*) 

o 

(P 2 , «V 2 -) * (Q**, 6fi a + ) ♦ (P.*, 6?,+) - (02, «E 2 ') 
where for the convenience of notation, tie shall use interchangeably the notation 
(X,T) 5 xTy 

Using these relations in equation (15), we obtain 

6J * - (Q lf 6R-j + ) - (P t , 6V!") ♦ (Q 2 , 6R 2 ") + (P 2 , fi? 2 +) 

♦ (On « ) - (Qb + , fiRg + ) 

The terns involving the vector P m in equation (15) have beat cancelled since the 
primer vector is continuous at the intermediate impulse; i.e., P B ~ = p B + = f b . 

To the first order, the variation in the position and velo£.ty vectors are 

6R = dR - Vdt 

$ = dV - Vdt 



( 16 ) 


Evaluating these equations at the times pp|+, t^ - , and t 2 ~ andfgisigg 

the results in equation (16), we obtain 


6J = - (Q-j^dRi + Pi^dV-j ) + (Q 2 ^dR 2 + P 2 ^dV 2 ) 

♦ (Q 1 ^V 1 + + P 1 T V 1 “)dt 1 - (Q 2 Tv 2 - + P 2 T V 2 +)dt 2 (17) 

+ - H m “ ) dt m - (Qni + - Qn“)^ dR n 


where H is the Hamiltonian function defined by equation (11). Equation (17) 
will be known as the generaSrdifferential cost function for the C/W equations. 

The first two terms on the right of this equation represent the variations in 
the cost function due to changes in the initial and final state, respectively. 

We shall have more to say about these two terms vhen we deal with specif i<^>ound- 
ary conditions. The third and fourth terms in equation (17) are the variations 
in the cost function because of changes in the initial and final times, respec- 
tively. If we are departing the initial orbit at the optimal time and arriving 
on the final orbit at the optimal time, and there are no other constraints, the 
coefficients of dt -j and dt 2 (respectively) will be zero. The last two terms 
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in equation (17) are the variations in the cost function because of changes in 
the intermediate impulse times and position vectors, respectively. For this lin- 
ear problem, there can be no more than four intermediate impulses (ref. 10). 

Note that on an optimal trajectory, the Hamiltonian and the vector Q must be 
continuous across the intermediate impulses. The vector P has already been 
required to be continuous by virtue of its definition (eq. (9)). 

As previously stated, equation (17) defines the general differential cost func- 
tion for the C/W equations. From this expression, it can be determined how any 
two-impulse trajectory or any two-impulse segment of an N- impulse trajectory can 
be improved: by changes in (1) the initial state vector or time, (2) the final 

state vector or time, or (3) the intermediate impulse times and position vectors. 


4.1 CRITERIA FOR AN INTERMEDIATE IMPULSE 

Consider a reference trajectory J, consisting of two impulses between fixed 
boundary conditions. Consider a perturbed trajectory J' between the same 
boundary conditions but consisting of three impulses - the intermediate impulse 
occurs at a time tj, and at a position vector R m + 6R m . The vector 6R m is 
the perturbation from the reference trajectory at the time tg,. We shall 
attempt to use a comparison between these two trajectories to determine under 
what conditions will the perturbed trajectory (the one with three impulses) 
have a lower cost than the reference trajectory. From equation (14), the 
difference in cost between the two solutions can be expressed as 


dJ = J' - J = 6J = |6v m + - 6v m -| - P 2 t 6V 2 - (18) 


since 6V-|“ and 6V 2 + are zero by definition of the boundary conditions. 
Evaluating the adjoint equation (eq. (13)) on the two segments of the perturbed 
trajectory, we have 


(P 1f «V! + ) = (Qnf, 6R m -) + (P m “, 6V m -) 
(P 2 , «V 2 -) = (V. «V> + (Pm + » «V n + ) 


since = 6R 2 =0 by virtue of the boundary conditions. 

Using these relationships in equation ( 1 8 ) , we have 

5J s (Qm i 6R m ) ~ (Qm + » 6hju‘> + |^m + “ i 
♦ (P m -, 5V m -) - (P ra + , 5V) 


(19) 


10 
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But the oostate vectors P and Q are evaluated on the reference trajectory 
and hence are continuous, or 


p nT s R m + s p m 

Qm~ = Qm'*’ 5 Qm 


Also, the variations In the position vector at the Intermediate times are 


s dR m - Vjjj _ dtjjj 


= dR m ” Vm^dtjj 

Using these results in equation (19), we obtain 

6J = |<5V m + - 6V m -| - P m T (6V m + - 6V m -) (20) 

since dt m is zero by definition of the perturbed trajectory. 

Equation (20) is homogenous in the intermediate impulse vector 6V m + - 6V m " 
for if we define a scalar u and a unit vector L as 


v = |<5V - 6v n -| 


L 


«V m + - «V m - 


v 


then equation (20) may be expressed as 


6j = v (1 - P m T L ) (21) 

The criteria for an intermediate impulse can now be established: if |P m | > 1, 
then 5j < 0, and J > J* ; or, the reference trajectory cost is greater than 
the perturbed trajectory cost. The reference trajectory cost can be improved 
by applying an intermediate impulse at the time t^ in the direction of P m . 
The greatest decrease in cost will occur if the intermediate impulse is applied 
at the time when |P m | is a maximum. 


11 
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Finally » if a reference trajectory exists for which |P_| > 1 at some tiae t^, 
how much should this trajectory be perturbed, and in whioh direction, suoh that 
the perturbed trajeotory has a lower oost? From equation (3). evaluated on the 
two segments of the perturbed trajectory, we have 

S^(t 1 ) a A(ti, t*) S-(t n ) 

S-(t 2 ) * A(t 2 , t m ) S+(t m ) 

Partitioning the matrix A as 

♦ ll ^12 
A a 

$21 $22 

The position vectors on the two segments can be expressed as 

R 1 = <hl(tl» W R m + $12 » W v m“ 

r 2 = ^*11 C 1 2 * tjj) R n» + 4* 12 (^2» ^m) ^m + 

Solving for the vector V [: + and V n ", we have 

V = 4> l2 - 1 (t 2 , t m )(R 2 - <t>n(t 2 , W R m ) 

V m ” = ^12”^^1» ^m^ R 1 “ 4> 1 1 ( 1 1 , t m ) R a ) 

where it can be ascertained that 4> 12 ” ^ exists if ujt t n2TT, n a 0,1,... 
Subtracting the two velocity vectors at the time t m , we obtain 

Pm 

V m + - V m " a AV m = p 2 - Pi ♦ MR m a V — (22) 


where 
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P 1 * ta) Ri 

p 2 s ^m) r 2 

and the matrix M is 

M s ♦i2“ 1 ^ t 1* tm) ♦llCti, ta) - <t>i2 _1 (t 2 , *m) 4>ll(t 2 , tm ) 


Solving for the intermediate position vector R m from equation (22), we have 


M-1 


m 


ra 


♦Pi - P2 


(23) 


The only unknown is this equation in the scalar v - the magnitude of the inter- 
mediate impulse. A second-order technique could be used to obtain an accurate 
value for V; however, experience has indicated that the actual value is not sig- 
nificant to the final solution, as long as it remains small relative to the 
cost J. 


4.2 BOUNDARY CONDITIONS 


4.2.1 Fixed-End-Conditions, Two-Impulse Transfer 

Consider the problem of finding the optimal two-impulse transfer trajectory be- 
tween fixed initial and final state vectors. Because intermediate impulses are 
excluded, t B and R m are fixed and therefore, 

dt m = dR m s 0 

Also, because the initial state vector (R i , Vi") and the final state vector 
(R 2 , v 2 + ) are r’ixed, we have 

dRi = dVi = dR 2 = dV 2 = 0 

Using these results in equation (17), the differential cost function reduces to 
6J = H^ti - H 2 dt 2 
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But beoause we are dealing with a two-impulse trajeotory, the Hamiltonian la con- 
tinuous, or 


H 1 = H 2 = H 


Also, if we define the transfer time as 


T = t£ - t-j 

the differential oost function for a time-open, two-impulse transfer between 
fixed initial and final states is 


6J s -Hdx 


(24) 


Equation (24) indicates that the oost function J will be an extremal when 
H s 0. Since H is defined in terms of elementary functions (eq. (11)), the op- 
timal time to transfer can be readily computed. The results of this solution 
will agree with those obtained in reference 11. 


4.2.2 Time-Open, Orblt-to-Orblt Transfer 

Consider the problem of transferring between two orbits, both fixed in shape and 
orientation, in which the transfer time and the departure and arrival times from 
the initial to the final orbits (respectively) are free. What are the necessary 
conditions for optimality to be satisfied for this solution? 

The differentials of the initial and final state veotors can be expressed as 


dRi = V'-dxo 
dV," = VrdTQ 


dR2 = V 2 +dxF 
dV 2 + = V 2 + dtp 


where to and Tp are time measurements in the initial and final orbits, re- 
spectively. Using these results in equation (17), the differential cost func- 
tion can be expressed as 


6J = -HrdTo + H 2 + dTp ♦ H^dt! - H 2 *dt 2 
+ (Hjn + “ Hjj, )dt m - (Q m + - Qjji )^ dR m 


(25) 
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Since all the differencials in this equation are independent and nonzero, 
the oost function J will be an extremal when all the ooeffioients of the 
differentials are zero. For an N-im pulse problem, there are 4(N-1) equations 
in terms of 4(N-1) unknowns. The optimal solution can be readily computed 
since all the ooeffioients are expressible in terms of elementary funotiona. 


4.2.3 Rendezvous 

As a final example, let us determine the neoessary conditions for an optimal, 

N- impulse, C/W rendezvous trajectory. This problem can be characterized 
as a time-open transfer; however, one that is funotlonally dependent on the 
motion in the initial and final orbits. That is, the ooasting times in the 
Initial and final orbits, to and tp (respectively) are not independent of the 
departure and arrival times, t-\ and t 2 (respectively), but are constrained by 
the relationships: 


dio = dti 


drp a dt2 


Using these relationships in equation (25) and the definition of the Hamiltonian 
from equation (11), we obtain 


6J = Q 1 T(V 1 + - Vr) dt! ♦ Q 2 t (V 2 + - V 2 “) dt 2 
* (V " »„-) dt m - (V - Qn-jTdRn 


But the vectors V-| + - V^“ a | AV t | P i , V 2 + - V 2 " = | AV 2 1 P 2 and the generalized 
differential cost function for a C/W rendezvous is 


6J a | A V t | ( P t - j ) dt -j ♦ |AV 2 |(P 2 T Q 2 )dt2 
♦ (H m + - Hjjj ) dtjj - (Q m + - Q m ") T dR m 


( 26 ) 


Note that for an optimal departure and arrival time from the initial and final 
orbit (respectively), the primer veotor P and its companion vector Q must be 
orthogonal. The vector Q, however, is not the negative derivative of the 
vector P, as can be verified in equation (7). 


5.0 EXAMPLE PROBLEM 


As an example to illustrate the analysis and computational procedure, a sequence 
of rendezvous trajectories will be computed. The target T is assumed to be at 
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rest and at the origin (fig. 1) of a coordinate system looated in a oiroular 
orbit with an altitude of 267 n. ai. at a time t. The ohase vehiole C is 
at rest, (with respeot to the target) and displaoed 10 n. mi. along the negative 
y-axis at a time 0. Thus, the boundary oonditions for the problem are 


RT(0) s (0, -10 n. mi., 0) 
R t (t) * (0, 0, 0) 


VT(0) * (0, 0, 0) 
V T (x) s (0, 0, 0) 


Let us now examine a sequence of rendezvous trajectories between these boui.dary 
conditions for various rendezvous times t. In figure 2, the cost of these solu- 
tions is indicated by the ourve labeled ^(NC) (the sum of the magnitudes of 
the two impulses with no initial ooast). The costate variables P and Q, and 
specifically the primer magnitude |P|, for any one of these solutions can be 
computed from equation (7) and may appear as illustrated in figure 3. Prom equa- 
tion (9) i the magnitude is unity at eaoh of the impulse times and in general 
will be less than unity on the interval (0, t). However, if the coefficient of 
dt-| (eq. (26)) is oomputed, it will be determined to be nonzero and for this 
problem it will also be positive. This implies that the oost of the solution 
can be reduced by a negative coast of duration to computed from the zero of a 
function f given as 


f ( tq) = |AV n | (P-,T Ql ) - o 


In figure 2, the two-impulse cost function with an optimal initial cost JgOfC) 
is also plotted. Note that the cost has been significantly reduoed by the intro- 
duction of an optimal initial coast. In figure 4, the optimal initial coast 
time tq is plotted against the rendezvous time t. Note that the rendezvous 
time has been maintained; i.e., To ♦ t-| s t where t<| is the transfer time be- 
tween the two Impulses. 

The question now arises, is it possible to further improve the optimal two- 
impulse trajectory? If we examine the primer magnitude time history for one of 
these optimal two-impulse solutions; i.e., t = 1000 seoonds, the primer magni- 
tude will appear as illustrated in figure 5. The primer magnitude is unity at 
the two- impulse times t-| s 0, t£ * Ti = 1450.3 seconds and is noted to be 
greater than unity on the interval indicating the oriteria specifying (eq. (21)) 
that this solution oan be improved by adding an impulse at the time where the 
primer magnitude is a maximum (t m = 926.3 seoonds). The starting position 
veotor (or this Impulse is given in equation (23). Henoe, a three- impulse solu- 
tion oan be oomputed and the optimal values for the coast time, arc times, and 
intermediate position vector can be computed from the zeroes of the coefficients 
of the respective variables in equation (26). 

The optimal three-impulse trajectory cost with an optimal initial coast J 3 (WC) 
is also plotted (dashed curve) in figure 2. Note that the optimal three-impulse 
cost is always less than or equal to the optimal two- impulse oost, but that 
these solutions may not always exist. In figure 2, note that for rendezvous 
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times prior to 655 seconds , three-impulse solutions do not exist. For this prob- 
lem, when three- Impulse solutions did exist, the oost was s oonstsnt value of 
134.7 fps. The aro tines between the three Impulses are illustrated in figure 6. 
Note onoe again that the rendezvous time is maintained; l.e., T o ♦ T 1 * T 2 * T * 

Finally, a significant observation about rendezvous trajectories is that the 
large variations in oost with rendezvous time oan be essentially eliminated by 
the introduction of optimal ooast and the optimal number of impulses. 
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Figure 1.- Positions of chase and target vehicles. 
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Impulsive arc times, t, and t 9 , sec 



Rendezvous time, r, sec 
Figure 6.- Optimal three-impulse arc time. 
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